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Abstract 

We offer a new form for the Sni term in the Hasselmann kinetic 
equation for squared wave amplitudes of wind-driven gravity wave. 

This form of Sni makes possible to rewrite in differential form the con¬ 
servation laws for energy, momentum, and wave action, and introduce 
their fluxes by a natural way. We show that the stationary kinetic 
equation has a family of exact Kolmogorov-type solutions governed by 
the fluxes of motion constants: wave action, energy, and momentum. 

The simple ’’local” model for Sni term that is equivalent to the ’’dif¬ 
fusion approximation” is studied in details. In this case, Kolmogorov 
spectra are found in the explicit form. We show that a general so¬ 
lution of the stationary kinetic equation behind the spectral peak is 
described by the Kolmogorov-type solution with frequency-dependent 
fluxes. The domains of’’inverse cascade” and ’’direct cascade” can be 
separated by a natural way. The spectrum in the universal domain is 
close to uj~^. 

1 INTRODUCTION 

The phenomenon of wind-generated gravity waves on the sea surface is a very inter¬ 
esting subject not only for oceanographers, coastal engineers and naval architects. 
The ocean is a great natural laboratory, and many phenomena taking place there 
are interesting for a broad community of physicists. Gravity and capillary surface 
waves on deep water represent the most conspicuous natural example of nonlinear 
waves in a strongly dispersive media. The statistical theory of such waves is called a 
theory of weak turbulence; it is an important part of general classical physics. This 
theory has been developing for more than forty years, and its basic concepts are 
now understood very well. However experimental data supporting this theory are 
scarce. The situation here is opposite to that of strong hydrodynamic turbulence 
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in an incompressible fluid. In this case we have a lot of experimental data but the 
theory is poor and inconsistent. 

The theory of weak turbulence is naturally applicable to the description of 
surface gravity and capillary waves on deep and shallow water. It can also be applied 
to acoustic Alfven waves in hydrodynamics, to many types of waves in plasmas, to 
waves in liquid helium and spin waves in ferromagnetics, as well as to Rossby waves 
in atmosphere and to internal waves; however experimental data collected so far 
in these areas of physics is too poor to perform a convincing comparison with the 
theory. Only the oceanographers who collect data on wave spectra, measured in 
the ocean, in lakes and in wave tanks for almost half a century, have accumulated 
vast quantities of very valuable experimental materials which could and should 
be compared with the predictions of the theory. (A contribution of professor M. 
Donelan in this process is really seminal, see for instance [Donelan, 1985]) This is 
an ambitious task; a lot of work must be done to perform it, and only first results 
are obtained in this direction. It was shown that the fetch dependence of wave 
energy and wave frequency, obtained in basic fetch-limited experiments, can be 
naturally interpreted in terms of the weak turbulent theory [Zakharov, 2002]. All 
experimenters agree that the observed spectra of wind-driven sea waves just behind 
the spectral peak have a universal power like form close to a;“^. We will show that 
this dependence can be easily explained in terms of the weak-turbulent theory. 

Once more, weak turbulence is the statistical theory of nonlinear waves in 
dispersive media. A central point of this theory is the following: a wave ensemble 
is described by the kinetic equation for square wave amplitudes. This equation has 
different names; for instance, the Boltzmann equation or the Hasselmann equation. 
Also, this equation could be called the Peierls equation, because it is nothing but the 
classical limit of the quantum kinetic equation for phonons, derived by Peierls and 
others in the late twenties. In this article we will use the term ’’kinetic equation”. 

The kinetic equation for gravity waves was derived by K. Hasselmann in 1962- 
1963 [Hasselmann, 1962; Hasselmann 1963]; now this equation is accepted as a basic 
model for the description of wave spectra evolution by the majority of oceanogra¬ 
phers. The kinetic equation in a truncated form (known as the DIA, or Direct 
Interaction Approximation) is widely used in the third generation models of wave 
prediction. However the physical effects described by the kinetic equation need to 
be commented. 

The main function of the nonlinear interaction term Sni in the kinetic equation 
is a very intensive redistributor of energy, momentum and wave action along the 
spectrum. Due to S'„;, the direct cascades of energy and momentum as well as 
inverse cascades of wave action are formed. These processes govern the evolution of 
the spectral peak and play a central role in the formation of the universal powerlike 
spectrum behind the spectral peak. In the simplest idealized cases, these spectra 
are Kolmogorov-type weak turbulent spectra, corresponding to constant fluxes of 
energy and wave action. Strictly speaking, such spectra are realized in physical 
systems, where domains of forcing and damping are essentially separated in A"- 
space. In the wind-driven sea the source of wave action is concentrated near the 
spectral peak, while the source of energy is distributed along the spectrum. On the 
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first glance this fact is an impediment for application of Kolmogorov-type theory; 
actnally this is not a serious difficulty. In a realistic situation the fluxes of energy 
and wave action are functions of frequency. However this does not affect essentially 
the spectral shape; it remains close to 

One more point is important. If the wind-driven sea is well-developed, then 
the main part of momentum fluxes from the wind is concentrated in short waves. 
This fact, experimentally established and essentially stressed by M. Donelan, can 
be naturally explained in terms of the weak turbulent theory. 

There is no reliable parametrization for the white capping dissipation. However 
it seems very probable that this fundamental process generates dissipation only 
in the short-wave region. Indeed, wave breaking makes the surface more smooth, 
acting like an efficient viscosity or even super-viscosity. Moreover, the wave breaking 
generates turbulence in the boundary layer with a thickness comparable to the 
length of the most dissipative waves. It is known that this layer is much thinner than 
the wave length of the spectral peak. Hence it is possible to suggest that the white 
cap dissipation in the area of the spectral peak can be neglected for the developed 
sea. In the first approximation the evolution of the spectral peak is described by 
the ’’conservative” kinetic equation such that includes only the time derivative, the 
advective term, and S„i. In the range of high frequencies the influence of wind 
forcing and white capping can be taken into account as a ” boundary condition”. 
This condition defines the flux of wave action to the long wave region. 

The ’’conservative” kinetic equation has a family of self-similar solutions de¬ 
pending on two free parameters. It was shown recently that by choosing the pa¬ 
rameters in a proper way one can explain the results of major fetch-limited exper¬ 
iments made during the last three decades, including the JONSWAP experiment 
[Zakharov, 2002]. We will present a detailed description of this study, supported 
by a massive numerical experiment, in an another article. In this paper we present 
the basic ideas of the weak turbulent theory using the simplest theoretical model 
of Sni known as ’’local” or ’’diffusive” approximation. 


2 Basic Theory 

We assume that the flow in the wave motion is potential v = V*!". The condition 
of incompressibility imposes on potential $ the Laplace equation: 

A$ = 0. (1) 

If T) is the shape of the surface, then equation (1) should be solved in the domain 
z < Tj under the boundary condition 


dz 


n = 

^ 0 . 
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In the linear approximation we should solve equation (1) in the half-space z < 0. 
The shape of the surface and the potential on the surface, ry and ip, are canonically 
conjugated variables; then the Euler equation for the potential flow of an ideal fluid 
with a free surface can be written 


Vt = 




I 


( 2 ) 


The solution of linearized motion equation (2) is the propagating monochro¬ 
matic wave 


/ 2u}r 

T] = W —^ Aq cosffcr - ujkt - (p), 

V 9 


2ujj: 

Ip = , Aosin{kf-uJkt - (p), 

V 1^1 


where uj^ = 


!g\^ and k is the wave vector. We can call 

A = Ao 


the complex amplitude of the propagating wave. The normalization of A is taken 
in such way that the energy density is 

E = wAp. 

By definition, Aq = E/oj is the density of ’’wave action”; then the sea surface is a 
composition of propagating waves 


9= [ \l^ (Ak + A*_^) dk, 

A real sea should be described statistically; to do this let us introduce the 
spectral density of wave action, assuming that 

{AuAl,) = gNj: 5{k-k). 

Then we can express the spatial correlation function 

F{R) = (’ 7 (^)’ 7 (^+^)) 
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in the form 


F{R) = j u)t Nt cos kRdk. 


In this case, the mean squared derivation a is given by the formula 




dk. 


Further, let us denote 7V^; this is the energy density in iti-space divided 

by g. It has dimension . Now we can express r]{r) through its Fourier transform 

dk 

vir) = J dk 

and define the spatial spectrum as follows 

{dkdk') = h 5{k + k), 

F{R)= [ I^e-^^Uk. 


Comparing with (3) we obtain 




a= Ij: dk. 


Further, it is convenient to introduce complex amplitudes 

H = ‘^^^k 

and derive the motion equation (2) in the form 

dat . SH 

k 

where the Hamiltonian iJ can be expanded in powers of 
H = H 0 + H 1 +H 2 + 


^0 = / Wr or dk, 


H, = 


2 J ^kkik 2 (®fc®fei®fc2 


* * 


x5{k — ki — ^ 2 ) dk dki dk 2 


3 J ^kkik2 


0^2 + ttr Or a? 

k ki k2j 


xS{k + fci + fe) dkdki dfe- 
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The Hamiltonian H 2 contains terms quartic in ai, a^. 

It is not very convenient to use equation (4). The cubic Hamiltonian Hi leads 
to the formation of ’’slave” waves; wave numbers and frequencies of ’’slave waves” 
are not connected by the dispersion relation. To separate ” slave” and ” free” waves 
one should perform a canonical transformation to new variables b^, eliminating the 
cubic term Hi. In new variables the Hamiltonian takes the form [Zakharov, 1999]: 

H = H 0 + H 2 , 

Ho = J ujj:bj:bldk, 

-^2 = ^ y^fcfci Safes ^S^Si^fes^fcs 

xS{k + ki — k2 — ko) dkdki dfe dfcs, 


and the dynamic equation 


% 

dt 


k 


in new variables is known as ’’Zakharov equation” [Zakharov, 1968] 


dbj: i 


Tu tt br bt bt 
fefcifcafes fci fe2 fes 


x5(k + fci — — ko) dki dk2 dko = 0. 


( 5 ) 


Explicit expressions for the coefficients of the Hamiltonian, the coupling coefficient 
^fcfciSafes’ canonical transformation can be found in (Zakharov, 1999). 

Equation (5), being approximate, has a very important feature: it conserves 
the total wave action, i.e., is adiabatic invariant 

iV = y dk, 


thus the kinetic wave equation is imposed to the correlation function of &-variables: 

{hbl>) = nj;Sik-k). 

On deep water we can put approximately 

nji-An^gNj:. 

It is important to stress that the kinetic equation describes not the real sea 
studied by experimenters but an idealized object: the ensemble of ’’free” waves 
filtered from the slave harmonics. On deep water the difference between these two 
ensembles is small (1 — 2%), while on shallow water the difference can be much 
more essential. 
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The kinetic equation in terms of N reads 


dN 


o 

+ ^+\/N = Snl+ S,n + Sds- 
Sk 


Here Sin and Sds are income from wind and dissipation due to white capping, and 
Snl has the form 


Snl 


IS {k, h,k,,h) 

- S[k + h-k2- h) 

xS (^LOj: + ujj:^ - dki dfe dks, 


where 


S (^k, ki,k2, fcs) = (27r)'‘7r5f^ 


kkik2 fes 
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can be found in [Hasselman, 1963], [Webb, 1978], [Zakharov, 1999]. 

The explicit expression for S is pretty complicated. The most important fact 
is the following: S (k,ki,k 2 ,k 3 ,'\ is a homogeneous function of sixth order 


S ( ek, eki, efe, efca, ) = e°S 


(k, ki, 


k2,k3 


For a rough estimate we can put 


S' ~ fc® ~ w 


12 


This is very fast growing function in frequencies. This fact is of a crucial importance. 
Most authors agree that Sin can be presented in the form 


Sin = l3{oj,9) N{k), 


where 


I3{uj,e) = ^ = 


to cos 9 
UJo 


Wo = 


UlO 


Here uio is the wind velocity at 10 meters height; = 0.1 ~ 0.3; Pa/Pw — 10“®; pa 
and pu, are densities of air and water. 

There is no agreement about the exact form of the function F{^). According 
to Hsiao and Shemdin [Hsiao, 1983] 


p = 0.12, F(0 = (0.85C-1)2, 


according to Donelan [Donelan, 1985] 


Aics 0.194, F(0 = (e-lf, 


while Tolman and Chalikov [Tolman, 1996] proposed a complicated form of F{^). 
In all these models /3(w) ~ w® as w —>■ oo. 
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According to Snyder [Snyder, 1981] 

M = 0.25, = 

In this case, (i{uj) ~ at large uj. 

An analytical expression for Sds is much less certain. Komen et al (1984) 
proposed the form 

= f—y (6) 

\apmj 

Here a is dimensionless steepness and Cdis is a dimensionless parameter. This 
formula is entirely speculative. It doesn’t have any theoretical foundation and is 
not derived from any real experiment in the ocean, lake or in a wave tank. Anyway, 
this formula is used widely in operational models (WAM, SWAM). It is supposed 
in most cases that 

n = 2, a = 3.33 x 10"^, apm = 4.5 x 10“^ 

In our opinion, expression (6) overestimates dissipation due to white capping in low 
frequencies. It can be used in absence of a better option; however on our opinion 
the parameter n should be essentially larger. If n > 3, the whole picture of ocean 
wave turbulence does not depend on a particular value of n. 


3 Constants of motion and their fluxes 


In this chapter we study the conservative homogeneous equation 


dN 


= S, 


nl 


(7) 


in absence of wind forcing and dissipation. It is considered that this equation has 
the following constants of motion - wave action, energy and momentum: 


N = J dk, 

E = J 

M = f kNj^dk. 


Formally speaking, this statement is correct; but the reality is much more compli¬ 
cated. Let us introduce polar coordinates |fc| and (j): 


\k\ 
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kdkd(j) = 

8 


—^ duj d(j), 
9 



and denote 


duj = Nk dk, 

2^3 

(j)) = - N 

9 


OJ 

9 


( 8 ) 


In what follows we understand N{uj, <j)) according to (8). We introduce also angle- 
independent spectra 


p2TT 

N{uj)= / N{uj,(j)) dcj), 

Jo 


E{uj) = uj 

= 


N{uj) 


2tt 


U) 
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N{u}, (j)) cos{<p) d<j). 


The conservative quantities in new variables take form 


^ = i 

= 


N{oj) doj, 
E{uj) duj, 

OO 

Mx{uj)duj, 


and conservation laws of these quantities can be written in the differential form 


dN{uj) 
dt ■ 
dE{uj) 
dt ~ 

dt 


dP 

dK 
doj ' 


Here Q is the flux of wave action to small wave numbers, while P and K are fluxes 
of energy and momentum directed to high wave numbers. A constant of motion 
is ’’real” if the corresponding flux is zero both at zero and infinite frequencies. 
Otherwise it is just a ’’formal” motion constant [Zakharov, Pushkarev, 2000]. 

Now let us introduce the differential operator 


_ 1 

^ 2^ 

and present kinetic equation (7) in the form 

dN{uj, (/>) 


1 d^ 
uj'^ d^P 


dt 


= LA. 


(9) 
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Here 


A{uj,(j)) = L ^ Snl, 

and H is a result of action on (jj) of a nonlinear integral operator 


A{u},(j)) = j F{u!,U!i,L 02,W3,4>- (j) 2 ,(l)- 4>i) 

xiV(wi, 4)1) N{u} 2 , 4>2) N{uJ 3 , 4)3,) 

xduJi duj2 duj3 d4)i d4)2 d4)3. ( 10 ) 

The explicit expression for F is given in Appendix. T’ is a homogeneous function 
of order 12: 

F{eu), ewi, ew 2 , eujs) = wi, W 2 , W 3 ) ~ 

Further, if we denote 

Aioo) = ^ J </>) d(/>, 

1 

B{u)) = — / B{uj,4>) cos 4) d4), 

271" Jo 

then the fluxes Q, P, K can be expressed in terms of A, B in the following form 



( 11 ) 

( 12 ) 

(13) 


Formulae (11-13) are of key importance for the theory of weak-turbulent spec¬ 
tra. 


4 Kolmogorov spectra 

In this chapter we study solutions of the stationary equation 

Snl = 0, 

which is equivalent to equation 

LA = 0. 

One class of solutions for (15) is given by solution of the equation 

A = 0, 


(14) 

(15) 
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and if these solutions exist, they are thermodynamic Rayley-Jeans spectra 

T 

'^k = -^—■ 

In the case of surface gravity waves these solutions do not exist because of the 
divergence of integrals in the operator A. To get physically significant solutions we 
can partially integrate equation (15) and put 

+ (16) 
CJ 

In this case, 


j4(cc) — loQ + -P, 

Biu;) = (17) 

UJ 


By substituting (17) into (11-13) we see that constants Q, P, and K in both cases 
are the same. 

Equation (16) defines the most general Kolmogorov-type solution of stationary 
kinetic equation; due to homogeneity of operator A this equation can be written in 
the form 




gsp-i 


R 




(18) 


with the energy spectrum 


E{u},(j)) 


ga Ps 


R 


ujQ 2kg 
~P' 'loP’ 


(19) 


Let us study the most important special cases. If Q = 0, KT = 0, formulae (18), 
(19) give the Zakharov-Filonenko Kolmogorov spectrum of the direct cascade 







( 20 ) 


Here 

Cp = i?(0,0,0) 

is the Kolmogorov constant of direct cascade (first Kolmogorov constant). We can 
offer another definition of Cp. 

Suppose that N{uj, (j)) is an isotropic powerlike function of w. 


N{uj)=uj-^. 


( 21 ) 
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Special consideration (which is not at home in this article) shows that integrals in 
A converge if 

19 

0 < a; < —. 

4 

Plugging (21) to (10) we obtain 

A(w) = (22) 


Apparently, 


fix—5 * 


Cf 


Spectrum (20) has a clear physical interpretation; this spectrum is a direct 
analog of the classical Kolmogorov spectrum of turbulence in an incompressible 
fluid. This spectrum is realized if there is a source of energy at small wave numbers 
and a sink of energy at high frequency region. 

The most general isotropic solution appears if AT = 0; then the spectrum is 




gipi 


(23) 


Function F{^) depends on one variable; obviously F{0) = Cp. If ^ oo, spectrum 
should not depend on P. Hence F{^) —>• CqjCp^!^ as ^ oo, equation (19) goes 
to Zakharov-Zaslavskii spectrum of inverse cascade 


A(w) 


11 1 
(jJ 3 


(24) 


where Cq is the Kolmogorov constant of the inverse cascade. Spectrum (23) pre¬ 
sumes that there is a source of wave action Q at high frequencies and sink of wave 
action at small frequencies. 

A general isotropic spectrum (23) is realized if there exists simultaneously a 
source of energy and sink of wave action at small frequencies together with energy 
sink and wave action source at high frequencies. 

If we study the most anisotropic case Q = 0, F = 0, then equation (16) has the 
following solution 




gt h{(j)) {Kg)l 

(jJ 3 


(25) 


From the symmetry consideration we can make a conclusion that 


h{(j)) = —h (tt — (/>). 


Hence solution (25) is not positive at all values. This is a reason to doubt that the 
general solution (19) is essentially positive and can be realized in the whole (to, </>) 
plane. Anyway, it can be used for approximation of real spectra in some finite part 
of wave-vector plane. 
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In the important case of Q = 0, solution (19) takes the form 




g-^p-^ 


H 



(26) 


Spectrum (26) can be found at small values of gK/uiP. Expanding in the Taylor 
series on this parameter, we obtain 




13 P 3 





a{4’)9K 

LoP 


The correction to the isotropic spectrum satisfies the linearized equation (16). As 
far as this situation is invariant with respect to rotations in the (w, </>) plane, 


a{(j)) = C 2 cos(j). 


Coefficient C 2 is known as the second Kolmogorov constant. 

If a; —>■ (X), then the right hand in (16) becomes independent of angle. This 
means that the Kolmogorov solution becomes isotropic at large uj. The real ob¬ 
served spectra remain anisotropic for arbitrary large oj. The explanation is the 
following: in real situation the momentum flux K is not constant but is approxi¬ 
mately proportional to frequency. 


5 Local diffusion approximation 


Many important features of the weak-turbulent theory can be understood in a 
framework of a very simple theory. 

Let us accept the following approximation for the operator A [Pushkarev, Za¬ 
kharov, 1999]: 

A{oj,4>) = ^N\ (27) 

9 

Here a is a dimensionless constant, which should be found by comparison with ex¬ 
periment. In this case the kinetic equation turns to the nonlinear diffusion equation 


dN a /1 52 1 92 \ 15 ^ . 

at ~ 54 (2 9^2 ^ w 2 9(),2 ) 


(28) 


Function A has a very simple form and a general Kolmogorov solution is 


Now 


E{uj,(j)) = 


g3 P3 


1-b 


w Q 2K g cos 


P 


ojP 


Cp — Cq — a 3 , 


(29) 
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Solution (29) is positive for all angles but in the case of large enough frequencies 
only; that satisfies the inequality 


LoP 

2kg 



> 1 . 


Comparing (26) with (9) we find that in this case 


9 


o 2 t: 


d<f), 


B{uj) = —OJ 
9 
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1 


('2'k 


cos <j)d(j), 


and for the general Kolmogorov solution we obtain 


A{uj,<j)) 

A{uj) 

B{u;) 


-P + w Q + 
P + uj Q, 

ui 


2Kg cos(j) 


Both Aiuj), B{uj) are essentially positive. This is correct for the general nonlocal 
case (12). The formula for A(u), (f) presumes that a solution has sources of energy 
and momentum, P and K, as well as a sink of wave action Q at the point a; = 0. 

In a framework of the local approximation we can efficiently study a forced 
stationary equation 

Snl + Sin + Sds = 0. (30) 

We can assume that 

Sin + Sds = /3(w, d) N{lo, 9), 

and restrict our consideration by the case of angular symmetry f3 = /3(w) only. 
Then equation (30) reads 




oj'-^ +l3{u})N = 0. 


2g^ duj^ 

Another form of this equation is the following: 


rA 


5^/3 / 3 (^) ^ 


p/3 ,.,5 

In a real situation solution N(uj) is concentrated in a finite frequency band 


(31) 


(32) 


iV > 0, LOi < UJ < (jJ2, 

N = 0, w < oji, UJ > UJ2- (33) 
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From continuity of N and dN/doj we obtain 


= 0 , 


= 0 , = 0 , 


Conditions (33) define the boundary value problem for equations (31), (32). 
Since in neighborhood of the ends of interval (33) there exists asymptotics 

A ~ -Pi(w —Wi)^ Pi > 0 
6 

A = -P 2 (w 2 —w)^ P 2 > 0 ( 

6 

we obtain from (32) the following expressions for /3: 


/3(cci) = -p2wf 


/3(W2)=-P2"w2^(^^j . (36) 

We can see now that a boundary problem has nontrivial solutions only if /3(w) is 
negative at both ends of interval wi < w < a; 2 . This conclusion is very general. 
To get a stationary solution of equation (30) we should have sinks both in low and 
high frequency regions. This statement without a proof can be found in the paper 
[Komen, Hasselmann, Hasselmann, 1984]. 

Condition (34) impose four restrictions on function N{oj) satisfying to a second 
order ODE. This is not too much because the ends of the interval uji < uj < 0 J 2 are 
unknown. They can be found from the following conditions for wave action and 
energy balance: 

PCJ2 

/ N{uj) duj = 0, 

J oil 

/ UJ f3{ui) N{uj) duj = 0. 

J uJi 

To satisfy the balance condition, we should have at least one domain of instability, 
where /3(w) > 0, inside the interval wi < w < uj 2 - In a typical situation there is one 
such area. In this case A(uj) has only one maximum at a point W 3 (wi < wa < UJ 2 ), 
and the whole interval could be divided into three domains: 

1 . Area, where A{uj) grows. Suppose, some interval A{uj) can be approximated 
by a linear function 

A{uj) = Q{uj - Wo). 

In this area, 

^ dA 

Q = = const, 


P = A-uj- 


—Wo <3 < 0. 
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This is the area of inverse cascade . A margin of this area is a frequency w*, where 
flux of energy P changes the sign P{uj*) = 0. 

2. Area near w ~ w*, where Aioj) is almost constant. Here Q is small, while 
P = A{uj*) is large and positive. This is the area of direct cascade . 

3. Area of dissipation, where dAjduj < 0. In this area P > 0, Q < 0. Both the 
energy and the wave action are carried out to a zone of high frequencies. 

In the area of direct cascade, equation (32) can be approximately integrated. 
We can rewrite this equation 


A 

dui 



+ 


g4/3 

Ai/3 a;4 


Ai/3 = 0, 


(37) 


and put 


dA 

w — < A. 

OOJ 


This makes possible to integrate (35); the integration yields 


In this approximation 


dA 

doj 


^ 9 ^ ^ 1/3 




2 5^/3 

3 ai/3 , 


/3(w) 


duj. 


P = .4 + P„l..(^) 


3/2 


Pn = 


d/2 \ 3^ 


3/2 


and Po is a slow function of w. For the spectrum in the direct cascade area we have 


P(w) 


1/2 


(38) 


Since in experiments w* < Wp (wp is a frequency of spectral peak), at the current 
level of experimental accuracy it is not easy to distinguish formula (36) from ZF 
spectrum 

We should stress again that the forced stationary equation (28) has a regular 
solution if and only if there are regions of intensive damping both in small and high 
wave numbers. What happens in other cases? Suppose, there is no damping at all. 
In other words. 


/3(cij) >0, Wl < W < W2, 

/3(u;) =0, w < wi., oj > UJ2 
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Then in the isotropic case, outside of the forcing areas the spectra turn to Kol¬ 
mogorov type spectra 

^ 4/3 pl/3 

= 1/3 4 ^ ^ > ^^ 2 , 

= ai/3cciV3 ’ 

and the fluxes of energy and wave action have the form 

pUJ2 

P= /3{ui) uj €{ui) doj, 

J U)i 

pUJ2 

Q= j3{uj) e{uj) duj. (40) 

J tO\ 

There is a difference of principle importance between the area of direct cascade 
w > W 2 and the area of inverse cascade w < wi. In the area of direct cascade, 
the integrals of motion, energy and wave action are finite. On the opposite, in the 
inverse cascade area the energy and the wave action diverge. We can say that the 
direct cascade has finite capacity, while inverse cascade has infinite capacity. 

A situation with no sink at high wave numbers is pure theoretical. Such a sink 
always exists due to pletora of physical reasons: viscosity, transformation of gravity 
waves to capillary waves, and finally due to wave breaking. A scrupulous consider¬ 
ation of these processes is not necessary for understanding: what happens near the 
spectral peak? Moreover, according to our preliminary study, Kolmogorov spec¬ 
trum of direct cascade can be formed in a finite time [Pushkarev, Resio, Zakharov, 
2000 ]. 

On the contrary, the Kolnogorov spectrum of inverse cascade, due to its infinite 
capacity, cannot be formed in a finite time. If there is no intensive enough damping 
at small wave numbers, the inverse cascade cannot be arrested. The downshift of 
spectral peak to small frequency area will continue infinitely until it will be stopped 
by topographical (better to say, geographical) factors. 

6 Discussion 

The first weak-turbulent Kolmogorov spectrum for gravity waves, ~ was 
derived analytically as an exact solution of the kinetic equation in 1966 [Zakharov, 
Filonenko, 1966]. The second Kolmogorov spectrum ~ was obtained 

in the same year in my PhD thesis in Novosibirsk [Zakharov, 1966]; in a regular 
journal the spectrum was published in 1982 [Zakharov, Zaslavskii, 1982]. It is called 
now Zakharov-Zaslavskii (ZZ) spectrum. 

For the first time, the spectrum ~ was observed experimentally by Toba 
in 1972. Since that this spectrum was observed by many researches [Forrestal, 1981; 
Kahma, 1981; Kawai et al, 1977; Donelan et al, 1985]. In 1987 Battjes et al found 
that spectrum fits the JONSWAP experiment much better than . In 1985 
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O. Phillips published a well known article, where he admitted that spectrum 
fits the experiment better than the ’’Phillips spectrum” However he did not 

offer a proper theoretical explanation of this fact. 

In 1982-83 Zakharov and Zaslavskii published in the Russian journal four arti¬ 
cles on application of the weak turbulence theory to the wind-driven sea [Zakharov, 
Zaslavskii. 1982; 1983]; soon after S. Kitaigorodskii used successfully these results 
for interpretation of experimental data [Kitaigorodskii, 1983]. Since that time the 
weak turbulent theory became known to the world community of oceanographers; 
however even now this theory is not completely accepted. For almost thirty years 
the obvious facts: 

1. Coj ~ uj~'^ is the exact solution of the stationary kinetic equation, 

2. ~ u)~'^ is the spectrum persistently observed in all experiments, 

coexist separately in the collective conscience, almost not interacting with each 

other. I believe, this is a unique situation in the history of science. 

The standard arguments against the weak-turbulent theory are the following 
(see, for instance, Komen, Cavaliery et al, 1994). 

1. Zakharov-Filonenko spectrum is isotropic, while the real spectra are anisotropic. 

This argument is not a serious one. The isotropic spectrum is the simplest 

example of weak-turbulent spectra. We showed in this paper that more general Kol¬ 
mogorov spectra, which carry momentum to high frequency region, are anisotropic. 
Anyway, they are very close to 

2. In the ’’classical” theory of turbulence for incompressible fluid, a source of 
energy is concentrated in small wave numbers while in the case of gravity waves 
the source of energy is distributed along the whole spectrum. 

One can add that the source of momentum is concentrated mostly in short 
waves. The answer is the following: in a real situation the flux of energy is not a 
constant inside the universal interval; this flux is a slowly growing function of fre¬ 
quency. This leads to nonessential modification: the appearance of slowly growing 
pre-factor proportional to (Inw/w*)^/^. 

We should stress again that the forced stationary kinetic equation (30) has a 
solution that contains a finite amount of energy if and only if there is an intensive 
dissipation in the low frequency region. This dissipation arrests the inverse cascade 
and is essential in the spectral peak area. Thus in this area equation (30) can be 
reduced to the form 

Snl + Sds = 0 . 

We like to stress that the physical origin of this low frequency dissipation is unclear; 
the very fact of existence of this dissipation and the whole concept of the ’’full- 
developed” are ’’mature” sea is questionable. We will discuss this subject in details 
in another article. 

In the universal region behind the spectral peak, the solutions of the full equa¬ 
tion (30) can be treated as solutions of the simple stationary equation 

Snl=0 

with frequency-dependent values of energy, wave action, and momentum flux. This 
is the central point of the theory of weak turbulence. We can add that this point is 
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supported now by massive numerical experiments (see for instance [Badulin et al, 
2002], [Lavrenov et al, 2002], [Pushkarev et al, to be published]). 

The research presented in this paper was supported by NSF grant NDM0072803 
and by the Army Corps of Engineers, RTDIE program, grant DACA 42-00-C0044. 


Appendix 

Solution of the equation 
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is given by the Green function 


A(uj, <p) 



G{u}, oj', (j) — (p') (j)') doj' dij)', 


(41) 


where 


Here 


and 


1 X A 's/iO iO 


n= —oo 


-1 


Ar,, = 


2n2 


1 ^ > 0 
0 ^<0 


0(0 = 

Equation (10) appears after substituting of S'„;(u;', </>') as f{ui', (j)') in formula (40). 
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